Jamming phase diagram for frictional particles 
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The non-equilibrium transition from a fluid-like state to a disordered solid-like state, known as 
the jamming transition, occurs in a wide variety of physical systems, such as colloidal suspensions 
and molecular fluids, when the temperature is lowered or the density increased. Shear stress, as 
temperature, favors the fluid-like state, and must be also considered to define the system 'jamming 
phase diagram' [I, @, Frictionless athermal systems[l|], for instance, can be described by the 

zero temperature plane of the jamming diagram in the temperature, density, stress space. Here we 
consider the jamming of athermal frictional svstems [sl. lol, [lol. [Til, [l^, Il3| such as granular materials, 
which are important to a number of applications from geophysics to industry. At constant volume 
and applied shear stress[3, 0]j ^6 show that while in absence of friction a system is either fluid-like 
or jammed, in the presence of friction a new region in the density shear-stress plane appears, where 
new dynamical regimes are found. In this region a system may slip, or even flow with a steady 
velocity for a long time in response to an applied stress, but then eventually jams. Jamming in 
non-thermal frictional systems is described here by a phase diagram in the density, shear-stress and 
friction space. 
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PACS numbers: 

Our analysis is based on Molecular Dynamics simula- 
tions of a suspension of soft-core spherical grains en- 
closed between two rough plates at constant volume frac- 
tion (/), and subject to a constant shear stress a. Periodic 
boundary conditions are used in the other directions (de- 
tail on the system and on the numerical method are given 
in the Supplementary Information). When the Coulomb 
friction coefRcient /i is set to zero, this model reduces 
to an assembly of frictionless particles, which, at zero 
applied stress tr = 0, jam at the random close packing 
volume fraction^ (j)j{a = 0) ~ 0.64. For 4> > (j)j{0) the 
shear and the bulk modulus grow as increases. When 
the system is subject to a small applied shear stress a, it 
flows when <j> is less than the jamming threshold 4>j{a), 
while it responds as an elastic solid when (jj > (j)j[a). 

At cr = 0, the presence of friction is known to mod- 
ify the location of the jamming point, which becomes 
also dependent on the preparation protocol@,[l3, 11, 
Here we focus on cr > 0, and show that friction gives 
rise to new dynamical regimes as illustrated in the phase 
diagram of Figure [TJi, where the flowing properties of a 
granular system with friction coefficient = 0.1 are sum- 
marized in the inverse density, 4'~^ i ^-nd shear stress, cr, 
plane. The initial state is prepared in such a way that 
no frictional contacts are present 0|, a condition which 
can be experimentally realized via high-frequency small 
amplitude vibrations [6[. At low density, in the 'Flow' 
regime, the system flows and reaches a stationary veloc- 
ity. For (j) larger than a threshold (jij^ = (l)j-^{a,^), the 
system enters the 'Flow & Jam' region. Here the sys- 
tem first flows with a stationary velocity (reached after a 
transient), but eventually enters by chance a microscopic 
configuration which is able to sustain the applied shear 
stress, and jams. The 'Flow & Jam' region is limited by 



a jamming line 4>J2 = '/'J2 (o": A*)- Above 4>j.-^ steady flow 
is never observed, and the system jams after a small slip. 
This 'Slip & Jam' region is limited by the line 4).j^{a,^) 
above which the system does not slip, but responds as a 
solid to an applied external stress. Examples of the sys- 
tem time-course in these different regions are illustrated 
in Fig. S2. 

How to define quantitatively the three lines? In the 
'Flow & Jam' region the system stops flowing after a 
jamming conflguration has been selected. The average 
time tjam the system flows before jamming is longer the 
smaller the volume fraction, and diverges approaching 
the 'Flow' regime, where it is actually inflnite. is de- 
fined as the volume fraction where tjam diverges coming 
from above (see Fig. [Jd). The line (pj.^ is defined by the 
divergence of the shear viscosity 77(0, cr, /i), defined as the 
ratio between shear stress cr and shear rate Vg/h, where h 
is the distance between the two plates, and Vs{(f), cr, /i) the 
shear velocity (soj. The viscosity increases as a power law 
with the volume fraction <f>, and diverges as <f> approaches 
i.e. 77 oc (0J2 — (f))"'^ (see Fig. [TJ;). The exponent 7 
appears not to depend on the shear stress, while it de- 
pends on the friction coefficient. The line (f>j^ marks the 
end of the 'Slip & Jam' region, as detailed in the Sup- 
plementary Information. We have defined the slip as the 
residual displacement of the top plate in a stress cycle, 
which allows to separate the slip of the top plate from its 
displacement due to the elastic deformation of the sys- 
tem on jamming. After preparing the system, we slowly 
increase the stress to its final value a, and then decrease 
it to zero. Below ^jg, this cycle is irreversible and the 
initial and final position of the top plate differ by the slip 
distance AL, while above (f),]^ the cycle is reversible and 
AL = 0. 
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FIG. 1: Flowing regimes of a granular system under 
shear stress, a In the "Flow" region the system flows with 
a steady velocity; in the "Flow & Jam" region the system 
first flows with a steady velocity, but jams after a time tjam; 
in the "Slip & Jam" region the system slips of a distance 
AL, never reaching a steady velocity, and then jams. In the 
"Jam" region the system deforms as a solid as soon as the 
shear stress is applied. The jamming time tjam, the viscosity 
77 and the slip distance AL have been fitted by power laws, 

tjam ~ - <^jj-'^, 77 ~ i^J, -^)-\ AL ^ (0J3 - 

for any given value of a and /.i. The volume fractions at 
which tjam diverges, ri diverges, and AL vanishes are shown, 
respectively, as squares (■), circles (•) and diamonds (♦). For 
cr = 2 10^'' and jj. = 0.1, b shows the power law divergence 
of tjam in — with a — 1.75 at 4)j^ = 0.622, and c that 
of the viscosity in — (j) with 7 — 0.75 at (jijj = 0.625. d 
Dependence of 4>j-^, (pj^ and 4>j^ on fi at a — 2 10^^. They 
conincide at = and decreases as /x increases. 

As discussed in the Supplementary Information, the 
line and c/jj^ of the the diagaram of Fig. [1] do not de- 
pend on the preparation protocol, for instance if friction 
is considered when inflating the particles in the prepa- 
ration of the initial state, as in several studies[10!, Illl |. 
This is because these lines are determined from extrap- 
olations in the flowing state where correlations with the 
initial state are by definition lost. The line c/jjg may de- 
pend on the protocol. The one we have determined is an 
upper bound with respect to those generated by all other 
protocols. 

<^Ji 1 4'J2 s-^d ^Js depend on the friction coefficient fi. 
Fig. [TJi illustrates their dependence on /x for <t = 2 10~^. 
The three lines coincide at /Lt = 0, if possible dynamical 
effects are neglected 0, but decreases with fi in such a 
way that (jjj^ {a, fi) > (pj^ (cr, /i) > 0,/^ (cr, ^) . In the limit 
of high friction, each line reaches a plateau. The depen- 
dence of on /i is very small, and only appears at high 
(j) or CT, due to the presence of a plastic response of the 
system (see Fig. S4). 

Such a dependence on friction leads to a jamming 
phase diagram for frictional particles characterized by 
three axis: the inverse density, the shear stress and the 
friction coefficient. In this phase diagram, schemati- 
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FIG. 2: Jamming phase diagram for frictional system. 

The jamming properties of frictional systems are illustrated in 
a diagram with axis the inverse density, the shear stress and 
the friction coefficient. At zero friction the jamming phase 
diagram is characterized by a "Flow" and by a "Jam" region, 
while in the presence of friction two new regions appear: the 
"Flow & Jam" region and the "Slip & Jam" region. 



cally shown in Fig. [2 the surfaces (j)j-^{a^ fi), (j)j^{a,^) 
and (/)j3 (cr, /i) enclose the regions of different flow prop- 
erties. The surfaces collapse in a line in the zero-friction 
plane, which coincides with the zero-temperature plane 
of the jamming phase diagram of frictionless particles. 
The phase diagram of Fig. [2] clarifies the intuitive ex- 
pectation that, when a frictional system jams after flow- 
ing, then it is possible to unjam it not only varying the 
density or the shear stress, but also by changing the fric- 
tion coefficient (which depends on humidity, temperature 
as well as on the presence of lubricants PJ, ll6|). The 
smallest value of the density at which we found jammed 
states is 4)j^{a ^ 0) ^ 0.585, close the smallest value 
reported in the litterture for jammed states in absence of 
gravity In the presence of gravity, looser states have 
been found [l3|. 

The changes in the structural properties of the system 
which occur crossing (/> , and 4>j^ are shown in Fig. [3] 
for /X = 0.1 and a = 2 10"'^. We discuss here the volume 
fraction dependence of the mean contact number Z , of 
the normal pressure on the confining walls P, and of the 
shear modulus G. The measure of G, which is discussed 
in the Supporting Information, is only possible because 
the system, despite having and anysotropic microstruc- 
ture, is not fragileQ and responds elastically and almost 
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FIG. 3: Mechanical properties Plot of the mean contact 
number Z (a), of the normal pressure on the shearing plate 
P (b), and of the shear modulus G (c) as a function of (f>, for 
a = 2 and n = 0.1. The vertical lines mark (f>j-^ (dashed), 
(f)j2 (dotted), and 4}Jr^ (plain). Circles are measures taken 
when the system flows, while diamonds are measure taken in 
jammed configurations. Open circles in the range i/ij^-^jj are 
measures taken in the flowing regime for t < tjam{<i>), before 
the system jams. 



the jamming properties of particulate systems. An open 
question is to explain the presence of the flow and jam 
region, where flowing frictional systems subject to a con- 
stant shear stress suddenly jam. In constant volume sys- 
tems, it is known that when the shear rate is fixed, large 
normal stresses fluctuations occur 2^, 21 1, as a con- 
sequence of a frustrated dialtancv|22|. When the shear 
stress is fixed, as in the case considered here, strong fluc- 
tuations arc not observed in the normal stresses, but in 
the shear rate. Wc speculate that at constant a a sys- 
tem jams in correspondence of a fluctuation so large that 
the shear velocity vanishes (see Supplementary Informa- 
tions). Therefore, the flow and jam phenomenology at 
constant shear stress can be seen as the counterpart of the 
large normal stress fluctuations observed at a constant 
shear rate[li, 0,(23. As jamming occurs when a perco- 
lating cluster of particles builds up, this argument may 
be related to the k~core percolation model introduced to 
describe jamming of frictionlcss particles[2^, or to mod- 
els which explicitcly take into account the constraint of 
mechanical equilibrium on each grain 24, 2^. However, 
it must be considcrd that while in absence of friction the 
percolating cluster can only emerge due to changes in the 
control parametrs, in the presence of friction such a clus- 
ter may spontaneously emerge while the system is flow- 
ing. The role of temperature in the jamming of frictional 
particles should be also investigated[2l]. Large colloidal 
particles, with a size smaller than roughly l/.tTO, are in 
fact at the same time small enough for temperature to 
influence their dynamics, and large enough to be charac- 
terized by frictional forces. Dense colloidal suspensions 
have actually already shown to behave as dense granular 
systems 
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isotropically to small external perturbations. In the flow- 
ing regime (circles) Z and P increase with (/>, while G is 
zero. In the jammed regime (diamonds) Z, P and G 
are roughly constant for < (f)j^ , while they increase as 
power laws for (p > (f)j.^ , where a continuous transition 
occurs. Measures taken in the flowing state in the range 
4>Ji~4>J2 (open circles) are taken for t < tjam{4>), before 
the system jams. These results are i n q ualitative agree- 
ment with recent experimental results [18| which have also 
show that, contrary to the frictionlcss case, Z and P do 
not vanish at the random close packing volume fraction in 
the presence of friction. Compared to previous numerical 
studies!^, [13, m conducted at ct = 0, our findings clarify 
that there is a whole volume fraction range where fric- 
tional granular systems may have the same mechanical 
properties. This volume fraction range can be identified 
with a constant Z line of the recently introduced Z-(j) 
diagram 

We have shown here that friction strongly influences 
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NUMERICAL MODEL 
Interaction between grains 

We have performed Molecular Dynamics simulations of a monodisperse system of particles of 
mass M and diameter D, based on a standard model for the grain-grain interaction, which is the 
linear spring-dashpot model. Two particles i and j, in positions and r^, with linear velocities 
and Vj, and angular velocities tUi and ujj, interact if in contact. The interaction force has a normal 
and a tangential component. 
The normal component is given by: 

where A:„ is the elastic modulus of the particles, = — rj, Sij = D — n^j = r,y/|ry |, 
■^riij — [(vi — Vj) • nijjriy . The effective mass is usually me// = MM/2M —1/2 (but see below 
for the interaction with particles of the confining boundaries). The parameter 7„ is fixed is such a 
way that the restitution coefficient is e = 0.88. 
The tangential component is given by: 

Fty = -ktMt.. -7tTOe//Vt.^-, 

where Uty is the elastic tangential displacement, and Vt^^ = Vy — Vn^^. Uty, set to zero at the 
beginning of a contact, measures the shear displacement during the lifetime of a contact. Its time 
evolution is fixed by vt; , w,; and Wj, as described in L. E. Silbert et al. Phys Rev. E, 64, 051302 
(2001). The presence of tangential forces implies the presence of torques, Ty = — l/2ry x Ft . The 
shear displacement is set to zero both when a contact finish {5ij < 0), and the Coulomb condition 
|Ftij I < I^Fty I is always enforced. Here jj, is the coefficient of static friction. 

As described in Sec. „ particles are enclosed between two rough plates. Each plate is made by a 
collection of particles that move as a rigid object. The bottom plate is fixed, and its particles are 
therefore considered to have an infinite mass. The top plate has a mass equal to the sum of the 
masses of its particles. The masses of the confining plates enter in the calculation of the effective 
mass in the interaction law. 

We use the value of the parameters of L. E. Silbert ct al. Phys Rev. E, 64, 051302 (2001): 
kn = 2 10^, kt/kn = 2/7, 7„ = 50, jt/jn ~ 0. Different values of the friction coefficient are 
investigated. Lengths, masses, times and stresses are measured in units of do ~ D, ttiq = M, 

to = ^/M/kn, CTO = kn/D. 



Preparation protocol 

The packings are constructed by first randomly placing the particles into the system with small 
radii, is such a way that no particles touch. Molecular Dynamics fictionless simulations are then 
performed by quickly inflating the particles radii in the presence of a small viscous damping force, 
until the radii reach their final value. Then, the system is allowed to relax until the kinetic energy 
vanishes. Friction is switched on after this procedure. 

Using this procedure, the volume fraction above which only jammed (finite pressure) states are 
generated is (j}rcp ^ 0.645 [H.P. Zhang and H.A. Makse, Phys. Rev. E 72, 011301 (2005)]. (jjrcp is 
approached exponentially fast as the system size increases (unpublished). 

Introducing friction after the preparation of the system allows the easy creation of dense packing 
of frictional systems. Experimentally, these high density states are generated via more complex 
procedures which allow for the continuous breaking of frictional contacts. Typical examples are 
vertical tapping, continuous high-frequency small amplitude vibrations, or thermal cycling. 
The influence of the preparation protocol is described in Sec^ 
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Computational details 

We solve the equations of motion of the system, mici = J2j + Fty and luji = J2j '''ij ^'^^ 
a velocity Verlet scheme, with an integration timestep 6t ~ 10^'^. The system reaches its steady 
state after a time of the order of T = 10^ in all regions of the phase diagram, but for the 'Flow & 
Jam' region, as long as cr > 2 10~^ (the minimum value we have considered). In the 'Flow & Jam' 
region, the system jams after a mean time which diverges on increasing the density, and therefore 
the steady state can be obtained only in a small volume fraction range. In this region, we have 
performed simulations lasting up to a time T = 5 10^. The simulations have been performed on 
a number of computer clusters. In 24/i, we simulate approximately a time 10'^, depending on the 
number of particles. We have performed simulations lasting up to 50 days. For each considered 
(j), (T, IJ, point, we have performed at least 10 different runs, starting from different initial conditions. 
In the "Flow & Jam region" , where we need statistics to properly measure the mean jamming 
time, we have performed 100 runs for each considered cr, yU point. 



INVESTIGATED SYSTEM 

We have investigated via Molecular Dynamics simulations a system of grains confined between 
two rough plates, as illustrated in Figure SI. The vertical distance between the plates is fixed, and 
a shear stress a along x is applied to the top plate. 




Figure SI: the investigated system. Grains are confined between two rough plates (red parti- 
cles) at a fixed vertical distance. A shear stress directed along x is applied to the top plate, while 
the bottom one is kept fixed. 



System size 

Particles are enclosed in a box of dimension Ix — ly = 16D, and Iz — SD. Periodic boundary 
conditions are used along x and y. The size of the vertical dimension 1^ is chosen to be comparbale 
to that of recent experiments [D. J. Pine, contribution to KITP Program on Granular Physics, 
2005, unpubhshed; J.-C. Tsai and J. P. GoUub, Phys. Rev. E 72, 051304 (2005); K. E. Daniels 
and R. P. Behringer, Phys. Rev. Lett. 94 168001 (2005)]. We have investigated the effect of the 
finite size of our system investigating system with Iz up to 64_D. The effect of the system size is 
described in Sec. „ 
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Volume fraction 



The volume fraction is equal to the volume occupied by the grains divided by the volume of 
the container, i.e. (p = Nvq/Vq, where Vq = Ixlylz is the volume of the system, N is the number of 
enclosed grains, and vq = I/GttD^ is the volume occupied by a single grain. 

Here, we have defined the volume fraction introducing a term which takes into account the effect 
of the rough plates protruding into the system. Due to the boundaries, the volume accessible to 
the grains is not Vq, but V = Vq — AV, where AV is an unknown corrective term. Since AV is 
much smaller than Vq, we have: 



In order to estimate AT^ we have determined the maximum number of grains Nmax which is 
possible to enclose into the system at zero applied stress, using the protocol described above to 
prepare the initial packing. We have then fixed (j){Nmax) = (f'rcp and determined AV inverting 




Eq. SI 



Eq. SI. 
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DYNAMICAL REGIMES 



Fig. S2 shows in the upper (lower) panel the time evolution of the position (velocity) in the 
Flow, Flow & Jam, Slip & Jam and Jam (in the insets) regions for cr = 2 lO^'^. The response of 
the system in the different regions can be summarized as follows: 
Flow: the system flows with a steady velocity reached after a transient. 

Flow & Jam: the system reaches a steady velocity after a transient. However, after flowing for 
sometime, it suddenly jams. 

Slip &: Jam: steady flow is never observed. The system jams after a small inelastic displacement 
of the top plate. 

Jam: the system responds as an clastic solid to the applied stress. 
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Figure S2: Dynamical regimes of a sheared granular system. 
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DEFINITION OF 



We have defined the shp distance AL((/i, cr) as the displacement of the top plate in a stress 
cycle. After preparing the system, we slowly increase the shear stress to its final value tr, and then 
decrease it to zero. Figure S3 (top panel) shows the displacement of the top plate position as a 
function of the shear stress for tr = 2 10"'^ and ^ = 0.1. Different curves refer to different values 
of the volume fraction, as shown. At small 0, the initial and final position of the top plate do not 
coincide: the residual displacement is our measure of the slip, AL, as exemplified in the figure for 
(j) = 0.642. The bottom panel shows that the slip distance decreases as a power law as the volume 
fraction increases: this allow to define (pj^ (a, fi) as the volume fraction where AL vanishes. 



a 

0.001 




AL(^ = 0.642) 




Figure S3: Left: displacement of the top plate in a stress cycle. The stress is first increased 
to its final value cr, and then decreased to zero. The residual displacement is our definition of the 
slip AL. Right: for a fixed value of a, the slip decreases on increasing the volume fraction, and 
vanishes at a volume fraction ^j,, which depends on a and fi. The straight line is a power law 
AL = a{(f> - 0J3)^ with a ~ 10^ b ~ 1.6. 

At very high values of the shear stress it is not possible to define (j),;^ this way, as one 
finds AL > even at very high volume fractions. The reason is that at high shear stress the 
system behaves plastically: one finds AL > not because the system slips, but because it deforms 
plastically in the stress cycle. However, when this is the case the dependence of AL on </> shows a 
clear crossover from a slip-dominated regime to a plastic-dominated regime, as shown in Fig. S4. 
This crossover allows to define (f>j^ as the infiection point of AL{(f>). 

The crossover from the elastic to the plastic regime is due to the increase of the number of 
contacts that are broken as the strain increase. At small cr, the strain of the system is small, and 
contacts do not break. At higher cr, the strain of the system is large, and contacts break. When 
a contact breaks, the tangential force between the grains is irreversibly destroyed, which is the 
microscopic origin of the plastic response of the system. 

lo' 

10° 
10'' 

<]1C 

10"' 

"'0.63 0.64 0.65 0.66 0.67 0.68 




Figure S4: At small a AL vanishes at high 0, and (j>j^ is defined as the volume fraction where 
AL vanishes as illustrated in Fig. S3. On the contrary at high cr, AL does not vanishes, but shows a 
crossover from a a slip-dominated regime to a plastic-dominated regime. In this case (f)j^ is defined 
as the inflection point of AL(0). 
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MEASURE OF THE SHEAR MODULUS G 



To measure the shear modulus G we have apphed to a system jammed under the action of a 
shear stress a a perturbmg shear stress 5a. The only non-zero components of Scr are 5axx and 5ayy, 
fixed such as Sul^ + ScFyy = 5a'^. The perturbing shear stress is therefore conveniently expressed 
in terms of its magnitude 5a and oi 6 ~ aicta.'ii[5ayy / Suxx) ■ The shear modulus G is defined as 
limber— >o 5a /e, where e is the shear strain induced by 5a. This definition is appropriate as for small 
5a {5a < 10~'^cr) the response of a jammed system is clastic (the strain is proportional to the 
stress) and to a good approximation isotropic (e only very weakly depend on the direction of 5a 
with respect to that of a, see below). Fig. S5 show the displacement 5{r) ~ {5x,5y) of the top 
plate position for different values of the volume fraction (</> varies from = 0.626 to = 0.663, 
the smallest (j) corresponding to the largest circle). Each curve is obtained by first applying a 
perturbing shear stress with {9 = 0), and then by increasing 9 from to 27r. The figure clarifies 
that systems jammed under shear are elastic, as each curve describes a close path. 

This result also clarifies that, even though the mechanical rigidity of a system jammed under 
shear stress originates from an underlying force network which is highly anisotropic, and which 
builds up only because of the presence of an applied shear stress, yet the system behaves as an 
elastic solid when a small perturbing shear stress is applied. It is therefore not 'fragile' as recently 
speculated [M.E. Gates, J.P. Wittmer, J. P. Bouchaud, P.Claudin, Phys. Rev. Lett. 81, 1841-1844 
(1998)]. Wc have checked that a fragile behaviour shows up in the response to larger perturbations. 




6x/D 

Figure S5: Response of jammed system to a small perturbing shear stress. 



The data of Fig.S5 suggests that the system behaves isotropically. To check whether this is 
actually the case, we have investigated the parameter 



m 



[5x^{9)+5y^i9)] 



1/2 



5r 



5r 



(1) 



where 5r = {[5x'^{9) + 5y^ 



^(6*) measures how close to the mean behavior the system is 



at each value of 9. Figure S6 shows that \S,{9)\ < 4%, suggesting that the response of the system 
is isotropic to a very good approximation. 




Figure S6: Anysotropy in the response of jammed system to a small perturbing shear stress. 
Different em'ves refer to different values of the volume fraction. 
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FINITE-SIZE ANALYSIS 

In this section, we discuss the robustness of the jamming phase diagram for frictional particles 
described in the main text to variations of the system size. We have kept fixed the size of the 
system in the transverse directions, Ix = ly = 16D, and varied the vertical size Iz- We compare 
the results for = 8, which arc the ones described in the manuscript, with results obtained with 
Iz — 16 and = 32, as obtained for cr = 2 10~^ and /i = 0.1. 

Finite size effects at (pj^ 

At the jamming line (jjj^, defined as the volume fraction at which the 'slip' vanishes, structural 
quantities have cusps, as shown in Fig. 3 of the manuscript. To investigate the dependence of the 
line (/>j3 on the system size, we have studied the size dependence of the location of the cusp in 
the pressure on the system size. As shown in Fig. S7, the cusp always occurs at the same volume 
fraction, implying that the line (jjj^ does not depend on the size of the system. This is not a 
surprise because (at small cr) ipj^ coincides with the random close packing volume fraction, which 
approaches exponentially fast its asymptotic value with the size of the system (when packings are 
generated with the protocol described in Sec. ^. 

0.01 
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0.006 

P 

0.004 
0.002 

(9.56 0.58 0.6 0.62 0.64 0.66 0.68 

Figure S7: Normal pressure acting on the top confining plate as a function of the volume 
fraction, for different system sizes. 
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Finite size effects at 

For each value of Iz, we have measured the shear viscosity 77 in the steady state, which appears to 
diverge as a power-law as the density increases. As shown in Fig. S8, data obtained with different 
sizes can be reasonably scaled on the same curve indicating that our system so large enough that 
finite size effects are small. The number of particles at 4).j^ varies between 2400 and 9600, depending 
on the size 

As (f) approaches , one enters the 'Flow & Jam' region of the phase diagram, where the system 
jams after flowing in a steady state for a time tjam- Since tjam becomes smaller and smaller as 
approaches 4>j^ , it is not possible to obtain reliable steady state shear viscosity data very close to 

(l>J2- 




2 

Figure S8: Log-log plot of the inverse shear viscosity 77"^ versus (/)j2 — (j), for different system 
sizes. The data collapse on the same master curve (77"^ ~ (^jj — 0)''', 7 ~ 0.85), indicating that 
finite-size effects are negligible. 
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Finite size effects at 

The jamming volume fraction <j)j-^ is that where the time ijam a system flows in a steady state 
before jamming diverges on decreasing the volume fraction. Since its definition involves a diverging 
time scale, its numerical identification is difficult, as well as the understanding of its dependence 
on the size of the system. To check for the presence of finite-size effects we have computed, for any 
given value of Iz, the probability p that a simulation jam in a given time T as a function of 0. The 
probability is computed over 100 runs which differ for the initial conditions, while the simulation 
time is fixed to T = 100. The results are shown in Fig. S9. 

1 



0.8 



0.6 

P 

0.4 
0.2 

8.62 ^0.622 '"' ' 0.624 0.626 

Figure S9: Fraction of simulation (over 100) that jams in a time T — 100 as a function of the 
volume fraction, for different system sizes. 

As the system size increases, at any given value of the fraction of runs which jam in a 
time T decreases. However, considering that we have only investigated a finite time T, and that 
the time required for a system to jam is expected to grow with the system size, one cannot 
draw from Fig. S9 any conclusion regarding the behavior of the line in the infinite system 
size, infinite T limit. Nevertheless, the flow & jam phenomenology appears to be relevant, as we 
observe it in systems with a size comparable to that of many granular experiments. We expect 
the phenolenology to be actually more apparent in experiments, where one can investigate a time 
T much larger than the one accessible in molecular dynamics simulations. 

Finally, we note that the flow & jam phenomenology is possibly related to the giant stress 
fluctuations observed in granular systems sheared at constant rate and constat volume fraction 
[see, for instance, B. Miller, C. O'Hern, and R. P. Behringcr, Phys. Rev. Lett. 77, 3110 (1996)]. 
At constant shear stress we do not observe large stress fluctuations, but large velocity fluctuations, 
the system jamming when there is a fluctuation so large that the velocity vanishes. In both cases 
the large fluctuations are observed in systems large enough to be of practical interest, while it is 
not clear if the fluctuations persist in the infinite system size limit. 
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INFLUENCES OF THE PREPARATION PROTOCOL 



Due to the presence of frictional forces, the response of granular systems to applied perturbations 
may depend on the particular protocol used to prepare the initial state. The phenomenology 
we have presented has been observed using a preparation protocol in which frictional forces are 
introduced after the system has reached a state of zero kinetic energy at the desired volume 
fraction. Our initial state is therefore memoryless. This protocol allows to access the whole zero 
pressure jamming phase diagram [C. Song, P. Wang and H.A. Makse, Nature 453, 629 (2008).] 
Experimentally, the expectation is that our initial states are those obtained compacting a granular 
system, for instance via high-frequency small- amplitude vibrations, which are able to destroy 
frictional contacts [G-J Gao, J. Blawzdziewicz, CS. O'Hern and M. Shattuck larXiv:0907.2106l 
(2009)], or via similar procedures. 

Here, we consider how our findings change when the initial packing is prepared using a different 
and popular protocol [H. P. Zhang and H. A. Makse, Phys. Rev. E 72, 011301 (2005); E. Somfai, 
M. van Hecke. W.G. EUenbroek, K. Shundyak, W. van Saarloos, Phys. Rev. E 75 020301 (R) 
(2007); and many others], where friction is always taken into account. Grains, initially placed in 
random positions with small radii, are inflated until they reach their final size. During inflation, 
frictional contacts are taken into account. We use the same inflation rate P both when using the 
protocol considered in the manuscript ('no friction protocol'), as well as when using the modified 
protocol ('friction protocol'). In Fig. SIO we compare, for cr = 2 10~^ and fj, = 0.1, the velocity of 
the shear plate (upper panel), and the pressure (lower panel) as obtained using the two protocols. 
The pressure is normal force acting on the top plate divided by its surface. 

The shear velocity is the same regardless of the initial protocol, in agreement with the expectation 
that the flowing systems don't remember their initial state. Accordingly, the line where the 
viscosity diverges (the velocity vanishes), is protocol- independent. The same is true for the line 
(not show), which is determined from the divergence of the jamming time, also measured when 
the system flows. 

The pressure, which is shown in the bottom panel, has a cusp at (pj^ (at small a). Figure 
SIO, therefore, clarifies that the line cfij^ depends on the preparation protocol. However, the line 
obtained with the 'no friction' protocol used in this work is a special one, and in this sense has 
to be preferred over the others, as it is an upper bound with respect to all possible lines obtained 
using different preparation protocols. 




Figure SIO: Location of the jamming transition lines using different protocol to prepare the 
initial state. 0,/^ and 4>J2 ^re protocol independent, while the line (jjJs depends on the protocol. 
Our estimate for is an upper bound for all possible estimations one can obtain using different 
protocols. 
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DYNAMICS OF JAMMING 

In the 'Flow & Jam' region of the jamming phase diagram, a system flows for a time tjam{4>), but 
then suddenly jams. This is a peculiar fluid to solid transition since it is not driven by changes in 
the control parameters - it is the result of a fluctuation which brings a flowing system in a jammed 
configuration able to support the applied stress. 

In Fig. Sll, we show the time evolution of several quantities as the system flows and then 
jams, at ~ 0.622, a = 2 10~^, and fi = 0.1. In panel (a), we plot the velocity of the top 
plate, which reaches a (fluctuating) steady state after a transient, but vanishes after a long time. 
Panel (b) shows the time evolution of the translational kinetic energy {Kt = J2i=i ~ 
of the rotational kinetic energy {Kr = J2iLi = 1/2^^?) > of the elastic energy due the normal 
{Un = = l/2fco|i5ijP) and to the tangential {Ut — J2i^j — l/2fct|uijp) interaction. The ratio 

between Un/Ut is roughly 100, as expected as the tangential force is bounded by the normal one, 
|ft| < iMfnl; so that Un/Ut — |ftp/|f„P ~ 1/m^- When the system jams, the elastic energies reach 
a plateaux, while the kinetic energies vanish. Panel (c) shows the evolution of the mean contact 
number, which increases on jamming. Finally, panel (d) shows the evolution of the normal pressure 
on the shearing plate. The pressure at each instant is equal to the force acting on the top plate 
divided by its surface (we are not averaging the signal over time) . The pressure fluctuates about a 
consant value: its fluctuations are much smaller than that found at constant volume and constant 
shear rate [B. Miller, C. O'Hern, and R. P. Bchringer, Phys. Rev. Lett. 77, 3110 (1996)]. 



a b 




Figure Sll: A typical run in the 'Flow & Jam' region of the jamming phase diagram. 
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The analysis of these quantities, as well as that of the fraction of the sliding contacts (a contact 
slides when the Coulomb condition is enforced), shown in Fig.S12, clarifies that in the jamming 
transition under shear the system quickly transients from a highly dynamic regime, in which the 
kinetic energy is serval orders of magnitude greater than the potential energy, and almost all 
contacts are sliding, to a jammed one, where the kinetic energy and the fraction of sliding contacts 
vanish. 




Figure S12: Time evolution of the fraction of sliding contacts. 



